The Two-Dimensional Square-Lattice S=l/2 Antiferromagnet Cu(pz)2 (0104)2 



O 



u 

e 

O 
o 



> 

On 
r- 
<N 
(N 

O 
O 



% 



N. Tsyrulini'2, F. Xiao^, A. Schneidewind'^'^ P. Link"^, H. M. R0nnow^, 

J. Gavilano^, C. P. Landee^, M. M. TurnbulF and M. Kenzelmann^'^ 

(1) Laboratory for Solid State Physics, ETH Zurich, CH-8093 Zurich, Switzerland 

(2) Laboratory for Neutron Scattering, ETH Zurich & Paul Scherrer Institute, CH-5232 Villigen, Switzerland 

(3) Department of Physics, Clark University, Worcester, Massachusetts 01610, USA 

(4) Forschungsneutronenquelle Heinz Meier- Leibnitz (FRM U), D-85H1 Garching, Germany 

(5) Institut filr Festkorperphysik, TU Dresden, D-01062 Dresden, Germany 

(6) Laboratory for Quantum Magnetism, Ecole Polytechnique Federale de Lausanne (EPFL), GH-1015 Lausanne, Switzerland 

(7) Garlson School of Ghemistry and Biochemistry, 

Glark University, Worcester, Massachusetts 01610, USA 

(8) Laboratory for Developments and Methods, Paul Scherrer Institute, GH-5232 Villigen, Switzerland 

We present an experimental study of the two-dimensional S=l/2 square- lattice antiferromag- 
net Cu(pz)2 (0104)2 (pz denotes pyrazine - O4H4N2) using specific heat measurements, neutron 
diffraction and cold-neutron spectroscopy The magnetic field dependence of the magnetic ordering 
temperature was determined from specific heat measurements for fields perpendicular and paral- 
lel to the square-lattice planes, showing identical field-temperature phase diagrams. This suggest 
that spin anisotropics in Ou(pz)2 (0104)2 are small. The ordered antiferromagnetic structure is a 
collinear arrangement with the magnetic moments along either the crystallographic b- or c-axis. The 
estimated ordered magnetic moment at zero field is mo = 0.47(5) /is and thus much smaller than 
the available single- ion magnetic moment. This is evidence for strong quantum fluctuations in the 
ordered magnetic phase of Ou(pz)2 (0104)2- Magnetic fields applied perpendicular to the square- 
lattice planes lead to an increase of the antiferromagnetically ordered moment to mo = 0.93(5) /xb 
at /ioH = 13.5 T - evidence that magnetic fields quench quantum fluctuations. Neutron spec- 
troscopy reveals the presence of a gapped spin excitations at the antiferromagnetic zone center, 
and it can be explained with a slightly anisotropic nearest neighbor exchange coupling described by 
J7 = 1.563(13) meV and J? = 0.9979(2)Jf . 

PACS numbers: 75.45. +j 75.30.Ds 78.70.Nx 



I. INTRODUCTION. 

Low dimensional quantum magnets are of great funda- 
mental interest. Unlike three-dimensional magnets, they 
support strong quantum fluctuations which can result in 
novel quantum excitations and novel ground states. Case 
in point is the antiferromagnetic S=l chain whose ground 
state features hidden quantum order and is separated by 
a finite energy from excited state^^^. In contrast, anti- 
ferromagnetic S=l/2 Heisenberg chains are gapless and 
feature fractionalized spin excitations as the hallmark of 
quantum criticality^. 

Increasing the dimensionality of a quantum magnet 
from one to two dimensions generally reduces effects 
of quantum fluctuations. The ground state of S=l/2 
square lattice Heisenberg antiferromagnet (AF) adopts 
Neel long-range order at zero temperature. Neverthe- 
less, strong quantum fluctuations arising from geometri- 
cal frustration may destroy long-range order in two di- 
mensions. 

Numerical studies of the two-dimensional (2D) S=l/2 
Heisenberg AF on a square lattice using quantum Monte 
Carlo, exact diagonalization, coupled cluster as well as 
series expansion calculations reveal a quantum renor- 
malization of the one-magnon energy in the entire Bril- 
louin zone and the existence of a magnetic continuum 
at higher energies^^. In recent years, quantum renor- 



malization effects at zero field have been studied using 
neutron scattering in a number of good realizations of 
S=l/2 square-lattice Heisenberg AF^™!^. The addition 
of antiferromagnetic next-nearest neighbor (NNN) inter- 
actions destabilizes the antiferromagnetic ground state 
and increases quantum fluctuations: according to the 
Ji — J2 modeP^H^, where Ji and J2 are the nearest 
neighbor (NN) and the NNN exchange interactions, re- 
spectively, different ground states are stabilized as a 
function of J2/J1. A possible spin-liquid phase appears 
to be the ground state for 0.38 < J2/J1 < 0.6 and 
collinear order was found for J2/J1 > 0.6. Our previ- 
ous stud}B^ of Cu(pz)2 (0104)2 has shown that even a 
small J2/J1 — 0.02 ratio enhances quantum fluctuations 
drastically, leading to a strong magnetic continuum at 
the antiferromagnetic zone boundary and the inversion 
of the zone boundary dispersion in magnetic fields. 

Here we present an experimental investigation of the 
2D organo- metallic AF Cu(pz)2 (0104)2, a good realiza- 
tion of the weakly frustrated J2/J1 — 0.02 quantum AF 
on a square lattice with Ji ~ 1.56 meV. Due to the small 
energy scale of the dominant exchange interaction, mag- 
netic fields available for macroscopic measurements and 
neutron scattering allow the experimental investigation 
of this interesting model system for magnetic fields up 
to about one third of the saturation field strength. We 
combine specific heat, neutron diffraction and neutron 
spectroscopy to determine the spin Hamiltonian and the 
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FIG. 1: (Color online) (a) Three-dimensional view of the crys- 
tal structure of Cu(pz)2 (€104)2 • The Cu^^ ions are shown 
as big spheres. The CIO4 tetrahedra are located between the 
copper layers and pyrazine molecules link Cu^^ in 6c-plane. 
The D atoms are not shown for simplicity, (b) The projec- 
tion of the crystal structure on the 6c-plane shows the Cu^^ 
square-lattice structure. The square lattice are shifted by 
(0,0.5,0) from one layer to the next. 
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FIG. 2: (Color online) Specific heat of Cu(pz)2 (0164)2 as a 
function of temperature for different magnetic field strengths. 
The measurements for magnetic fields parallel and perpen- 
dicular to the 6c-plane are shown in (a) and (b), respectively. 
For convenience the graphs are shown as lines connecting the 
data points and without error bars. 



key magnetic properties of this model material. Specific 
heat measurements show that the magnetic properties 
are nearly identical for fields applied parallel and perpen- 
dicular to the square-lattice plane. This shows that spin 
anisotropics are small in contrast to spatial anisotropics, 
and that it is sufficient to perform microscopic measure- 
ment for just one field direction. Our microscopic neu- 
tron measurements, on the other hand, provide informa- 
tion on the spin Hamiltonian that explain the nearly iden- 
tical HT phase diagrams for the two field directions. Spe- 
cific heat and neutron measurements of Cu(pz)2 (0104)2 
thus ideally complement each other. 

Deuterated copper pyrazine perchlorate 

Ou(pz)2 (0104)2 crystallizes in a monoclinic crystal 
structure described by space group 02/c, with lat- 
tice parameters a = 14.045(5)A, b = 9.759(3)A, 
c = 9.800(3)A and /3 = 96.491(4)°22. The crystal 
structure is shown in Fig. [I] The Ou^+ ions occupy 4e 
Wyckoff positions and pyrazine ligands link magnetic 
Ou^+ ions into square-lattice planes lying in the crys- 
tallographic 6c-plane. The Ou^+-Ou^+ NN distances in 
the 6c-plane are identical and equal to 6.92A^^. The 
two fold rotation axis (0, y, 1/2) and the mirror plane 
parallel to the ac-plane ensure that all NN exchange 
interactions between Ou^+ are identical. Tetrahedra of 
OIO4 located between the planes (Fig. flk) provide good 
spatial isolation of Ou^+ ions and substantially decrease 
the interlayer interactions. Thus, perfect square-lattices 
of copper ions with a superexchange path mediated by 
pyrazine molecules are formed in the 6c-plane (Fig. fib). 

The magnetic susceptibility shows good agreement 
with that of the 2D S=l/2 Heisenberg AF with an ex- 
change interaction strength of Ji = 1.53(3) meV. Small 
interlayer interactions result in a long range antiferro- 



magnetic order below Tn = 4.21(1) K. The ratio of 
interlayer exchange, Jx, to the dominant intralayer ex- 
change strength, Ji, was estimated as J±/Ji = 8.8 • 10~^ 

[22^241 



II. EXPERIMENTAL DETAILS. 

In order to obtain the field-temperature (HT) phase di- 
agram of Ou(pz)2 (0104)2 we measured the specific heat 
as a function of temperature for different magnetic field 
strengths using the Physical Property Measurement Sys- 
tem by Quantum Design. A single crystal of deuterated 
Ou(pz)2 (0104)2 with mass m = 13 mg was fixed on a 
sapphire chip calorimeter with Apiezon-N grease. The 
measurements were done using the relaxation technique 
which consists of the application of a heat pulse to a 
sample and the subsequent tracking the induced tem- 
perature change. The specific heat was obtained in the 
range from T = 2 K to T = 30 K in magnetic fields 
of up to jhqH = 9 T applied parallel and perpendicular 
to the copper square-lattice planes. The measurements 
were done with the steps of ATi = 0.05 K, AT2 = 0.2 K 
and AT3 = 1 K in the temperature ranges Ti = 2 — 6 K, 
T2 = 6 — 8 K and T3 = 8 — 30 K, respectively. Oare was 
taken to apply a small heat pulse of 0.1% of the tempera- 
ture step AT and each measurement was repeated three 
times to increase accuracy. Specific heat of Apiezon-N 
grease without Ou(pz)2 (0104)2 crystal was measured in 
the entire temperature range separately and subtracted 
as a background from the total specific heat of the sample 
and grease. 

The HT phase diagram and the ordered magnetic 
structure were studied by neutron diffraction using cold- 
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FIG. 3: (Color online) Neutron scattering peak intensity of 
Q = (0,1,0) as a function of temperature and magnetic field, 
applied perpendicular to the square-lattice planes. The re- 
sults obtained by specific heat measurements in magnetic field 
applied parallel and perpendicular to the Cu planes are shown 
by squares and circles, respectively. The neutron data, mea- 
sured for magnetic fields perpendicular to the copper planes, 
are shown by triangles. 



neutron three-axis spectrometer RITA2 at the Paul 
Scherrer Institute, Villigen, Switzerland. A crystal with 
dimensions 7 x 7 x 1.5 mm and mass of m = 85 mg 
was wrapped into aluminum foil, fixed with wires on a 
sample holder and aligned with its reciprocal [0, k, 1] 
plane in the horizontal scattering plane of the neutron 
spectrometer. Data were collected at T = 2.3 K and 
T = 10 K in magnetic fields up to /ioH = 13.5 T ap- 
plied nearly perpendicular to the [0, k, 1] plane using an 
Oxford cryomagnet. Measurements were performed with 
the pyrolytic graphite (PG) (002) Bragg reflection as a 
monochromator. A cooled Be filter was installed before 
the analyzer to suppress higher order neutron contami- 
nation for the final energy Ef = 5 meV. We also used an 
experimental setup without Be filter, which allowed to 
use the second order neutrons from the monochromator 
with Ei = 20 meV, thus allowing to access to reflections 
at high wave- vector transfers. 

The spin dynamics in the antiferromagnetically or- 
dered phase was measured using the cold-neutron three- 
axis spectrometer PANDA at FRM-2, Garching, Ger- 
many. Two single crystals with a total mass of m = 1 g 
were wrapped into aluminum foil, fixed on a sample 
holder with wires and co-aligned in an array with a fi- 
nal mosaic spread of 1°. Reciprocal [0, k, 1] plane of 
the sample was aligned with the horizontal scattering 
plane of the neutron spectrometer. These measurements 
were performed in zero magnetic field and at tempera- 
ture T = 1.42 K using a ^He cryostat generally referred 
to as an Orange cryostat. The final energy was either set 
to Ef = 4.66 meV or Ef = 2.81 meV using a PG(002) an- 
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FIG. 4: (Color online) (a) Neutron scattering intensity at 
Q = (0, 1,0) measured at T = 2.3 K and at T = 10 K as 
a function of a rotation of the sample around the vertical 
axis described by angle B. The inset (b) shows the neutron 
scattering observed at Q — (0, 0, 1) at the same temperature. 



alyzer. Data were collected using a PG(002) monochro- 
mator and cooled Be filter installed before the analyzer. 



III. RESULTS. 
A. Specific heat measurements. 

The temperature dependence of the specific heat of 
Cu(pz)2(C104)2 is shown in Fig. ^for different magnetic 
fields applied perpendicular and parallel to the copper 
square-lattice plane. At all fields, the temperature de- 
pendence of the specific heat reveals a well defined cusp- 
like peak, indicating a second order phase transition to- 
wards 3D long-range magnetic order. Previous zero-field 
studies of Cu(pz)2 (0104)2 did not show an anomaly in 
the specific heat^^. Most likely, the high accuracy of our 
measurements played a crucial role in detecting the zero- 
field anomaly in the specific- heat curve. The small size of 
the ordering anomaly is a consequence of the low dimen- 
sionality of the magnetism and an ordered magnetic mo- 
ment that, due to quantum fluctuations, is considerably 
smaller than the free-ion value. The HT phase diagram 
assembled from the specific heat measurements is shown 
in Fig. [3J The measurements show that the Neel tem- 
perature increases with increasing magnetic field, from 
Tn = 4.24(4) K at zero field to Tn = 5.59(3) K at 
/ioH = 9 T. 

We also observe an increase of the specific heat with 
increasing magnetic field in the paramagnetic phase just 
above the 3D ordering temperature. We propose that the 
field dependence of the specific heat data is a consequence 
of field-induced anisotropy in the 2D AF. In zero field, a 
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FIG. 5: The magnetic peak intensity of neutron scattering 
at Q = (0,1,0) as function of magnetic field measured at 
T = 2.3 K. The non-magnetic scattering was estimated from 
measurements at T = 10 K and subtracted from the overall 
peak intensity. The inset shows the ordered antiferromagnetic 
moment as a function of field. 



pure 2D Heisenberg AF orders at zero temperature, but 
quantum Monte Carlo simulations^^ have shown that ap- 
plication of an external field induces an Heisenberg-XY 
crossover and leads to a finite temperature Berezinskii- 
Kosterlitz-Thouless transition Tbkt^^^- One conse- 
quence of this crossover is the increase of Tbkt with ex- 
ternal field for up to i^ < J^sat/4 and then a gradual de- 
crease of the transition temperature with increased fields. 
While the zero-field 3D transition T/v in Cu(pz)2 (0104)2 
is driven by the combination of 3D interaction and in- 
trinsic XY anisotropy, the increase of T/v as a function 
of field may thus be driven by an increase of the effective 
anisotropy and the associated increase of Tbkt- Sim- 
ilarly, we propose that the increase of the specific heat 
above the 3D ordering temperature is caused by the field- 
induced XY anisotropy: In the 2D antiferromagnet on a 
square lattice, 2D topological spin- vortices appear above 
the Berezinskii-Kosterliz-Trousers (BKT) transition as 
the preferable thermodynamic configuration. In applied 
magnetic field the vortices unbind above the BKT transi- 
tion, leading to the increase of the specific heat above the 
ordering temperature. The anisotropy crossover thus af- 
fects the specific heat in a manner similar to the observed 
behavioi^^. 

Remarkably, the HT phase diagrams are identical for 
fields parallel and perpendicular to the square-lattice 
planes. This suggests that the dominant exchange inter- 
actions between nearest copper spins Ji in 6c-plane are 
close to the isotropic limit in spin space. This should not 
be confused with the strong spatial two-dimensionality 
ofCu(pz)2(C104)2. 



B. Magnetic order parameter. 

To determine the ordered magnetic structure of 
Cu(pz)2 (0104)2, several magnetic Bragg reflections were 
measured by neutron diffraction. Fig. |4] shows two mag- 
netic peaks measured above and below the transition 
temperature at Q = (0,1,0) and Q = (0,0,1) using a 
final energy Ef = 20 meV. This data directly demon- 
strates the presence of magnetic order below Tn- The 
magnetic Bragg peak widths are limited by the instru- 
mental resolution, confirming that the magnetic order 
is long-range. The field dependence of the magnetic 
scattering at Q = (0,1,0) measured using final energy 
Ef = 5 meV reveals an increase of magnetic intensity 
as a function of field from zero to /ioH = 13.5 T as is 
shown in Fig. [5j The magnetic scattering was determined 
by subtracting the non-magnetic background determined 
at T = lOK. The increase of magnetic diffraction inten- 
sity with field is most probably related to a quenching 
of quantum fluctuations by the magnetic field, that si- 
multaneously also leads to the observed increase of the 
transition temperature Tn- This result is in a good agree- 
ment with the specific heat data indicating enhanced XY 
anisotropy in the applied magnetic field. The intensity 
measured at Q = (0, 1,0) at T = 10 K as the function 
of applied field did not reveal any magnetic scattering, 
showing that magnetic fields do not lead to field-induced 
antiferromagnetic order in the paramagnetic phase. 

The critical magnetic behavior was studied by mea- 
suring the peak intensity of the neutron scattering at 
the antiferromagnetic wave vector Q = (0,1,0) and 
Q = (0, 3, 0) as function of temperature in magnetic field 
up to /ioH = 13.5 T. Typical scans are shown in Fig. [6] 
The solid line display that the increase of the antiferro- 
magnetic intensity in the ordered phase close to Tn is 
evidently steeper for high fields. The HT phase diagram 
compiled from the temperature scans is shown in Fig. [3] 
and it confirms the phase diagram obtained from specific 
heat measurements. 



C. Ordered magnetic structure. 

The symmetry of the ordered magnetic phase was stud- 
ied by neutron diffraction. Group theory was used to 
restrict the search only to magnetic structures that are 
allowed by symmetry. The magnetic Bragg peaks at 
Q = (0,1,0) and Q = (0,3,0) indicate that the mag- 
netic structure breaks the 0-centering of the chemical 
lattice and that Ou(pz)2 (0104)2 adopts an antiferromag- 
netic structure for T < Tn- Symmetry analysis revealed 
six basis vectors which belong to four irreducible rep- 
resentations and are listed in Tab. Ill (for details see 
Appendix A). 



The analysis is complicated by the fact that the single- 
crystal probably consists of two domains with inter- 
changed h- and c-axis, which are nearly identical in 
length. A twinning of the single-crystal in this manner is 
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FIG. 6: (Color online) (a) The temperature dependence of 
the neutron scattering peak intensity measured at the an- 
tiferromagnetic point Q =(0, 1, 0). The data collected at 
^oH = T, /ioH = 2 T /ioH = 6 T and /ioH = 13.5 T are 
shown by circles, squares, triangles and diamonds, respec- 
tively. The red lines are guides to the eye. The inset (b) 
represents the peak intensity at Q = (0, 3, 0) as the function 
of temperature obtained at zero field. Tn was found to be 
the same as for Q = (0, 1, 0). 



indicated by the observation of both the Q = (0, 2, 3) and 
Q = (0,3,2) nuclear Bragg peaks with similar intensity, 
although Q = (0,3,2) is not allowed for a C-centered 
lattice. 

The experimental data are consistent with both r2 
and r4 irreducible representations listed in Tab. [TT] and 



with two basis vectors 



and 



It is not possi- 



ble to distinguish between these two solutions because 



-> 



2 of the be crystallographic domain is identical to 
6 of the cb crystallographic domain, and the fits were 
made assuming a equal population of be and eb crystal- 
lographic domains. The ordered magnetic structure of 
Cu(pz)2 (0104)2 can have magnetic moments aligned an- 
tiferromagnetically either along crystallographic b- or e- 
axis as is shown in Fig. iTFa) and Fig. iT^b), respectively. 
Due to a small number of observed magnetic reflections 
and the crystallographic twinning, our experiment can- 
not distinguish between these two magnetic structures. 
The collinear spin arrangement in 6c-plane is consistent 
with the absence of the Dzyaloshinsky-Moriya interac- 
tions between NN. The spatial arrangement of the or- 
dered magnetic moments in adjacent square-lattice layers 
is ferro- and antiferromagnetic along ab- and ac-diagonal, 
respectively. This is consistent with the chemical struc- 
ture of Ou(pz)2 (0104)2, where the interlayer interaction 
pathway along ac-diagonal is shorter than the path along 
ab. 

The value of the ordered magnetic moment was ob- 
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FIG. 7: Two possible magnetic structures of Cu(pz) 2 (0164)2 
belonging to the second and the fourth irreducible representa- 
tions (see Tab.ro, are shown in (a) and (b), respectively. Two 
adjacent square-lattice Cu^^ layers, separated by a (0.5,0.5,0) 
lattice unit translation, are depicted by open and filled arrows. 
Cu^^ — Cu^^ interlayer interaction pathway along ac diagonal 
corresponds to the vertical distance between filled and open 
symbols in (a) and (b). 



tained from a minimization of ^ = IR, 



calc 



Rexpl, where 



Rexp is the measured ratio of the magnetic Bragg peak 
intensity to the nuclear Bragg peak intensity, Rcaic = 
|F(Q)magn|V|F(Q)nuciP, F(Q)magn and F(Q)nuci are the 
magnetic and nuclear structure factors, respectively. The 
fit was performed for two magnetic peaks observed at 
Q = (0,1,0) and Q = (0,3,0) and two nuclear peaks 
measured at Q = (0,2,4) and Q = (0,0,6). The ob- 
tained value of the ordered magnetic moment in zero field 
is mo = 0.47(5) jub- The comparison of Rcaic and Rexp for 
two magnetic and two nuclear Bragg peaks is presented in 



|i^(0,l,0)|^ 
|F(0,2,4)|2 


|F(0,1,0) ^ 
|F(0,0,6) 2 


|F(0,3,0) ^ 
|F(0,2,4) 2 


F(0,3,0) ^ 
F(0,0,6) 2 


Rexp X 10-^ 3.72(7) 


5.13(12) 


1.89(18) 


2.60(25) 


Rcaic X 10-^ 4.35 


4.57 


2.16 


2.26 



TABLE I: The measured and the calculated ratios of squared 
magnetic to nuclear structure factors for different Bragg 
peaks. The calculated values were obtained from a mini- 



mization of ^ = I Rcaic 
magnetic moment mo ■ 



Rexp I and correspond to the ordered 

:0.47/iB. 



the Tab. |Tj The calculated value of the ordered magnetic 
moment is smaller than the free- ion magnetic moment, 
indicating the presence of strong quantum fluctuations 
in the magnetic ground state of Ou(pz)2 (0104)2. The 
inset (b) in Fig. [s] displays the increase of the ordered 
antiferromagnetic moment from mo = 0.47(5)/iB in zero 
field to mo = 0.93(5) /iB in MoH = 13.5 T. This is direct 
evidence for the suppression of quantum fluctuations by 
the applied magnetic field due to induced XY anisotropy 
as suggested by our specific heat measurements. 
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FIG. 8: (Color online) A series of constant energy scans per- 
formed along the [0, k, 1] and [0, -1, 1] directions at different 
energy transfers AE in zero magnetic field and at T = 1.42 K. 
Please note the changing scale of the vertical axis for the dif- 
ferent scans. The solid lines correspond to a convolution of 
two Gaussians with the resolution function. 



D. Spin dynamics. 

The wave-vector dependence of the magnetic excita- 
tions has been measured using neutron spectroscopy. 
Constant energy scans were performed near the antiferro- 
magnetic zone centers Q = (0, 0, 1) and Q = (0, —1, 0) for 
energy transfer AE in the range from AE = 0.5 meV to 
AE = 3 meV and are shown in Fig.[8J The observed mag- 
netic peaks are resolution limited, indicating that these 
magnetic excitations are long-lived magnons associated 
with a long-range ordered magnetic structure. 

Constant wave- vector scans were performed at the an- 
tiferromagnetic zone centers in the energy transfer range 
from AE = meV to AE = 0.7 meV (Fig. [9^, b). These 
scans reveal a magnetic mode which is gapped and has 
a finite energy Ezc = 0.201(8) meV at the antiferro- 
magnetic zone center. The energy gap at the antifer- 
romagnetic zone center is attributed to the presence of 
a small XY anisotropy in the nearest-neighbor two-ion 
exchange interactions, because a single-ion anisotropy of 
type D(S^)2 is not allowed for S=l/2. 

Constant wave-vector scans away from the antiferro- 
magnetic zone center carried out at higher energies are 
shown in Fig. M^c). The energies of the magnetic exci- 



tation at the symmetrically identical antiferromagnetic 
zone boundary points Q^^i = (0? 0-5,1) and Q2b2 = 
(0, —0.5, 1) are equal to E^bi = 3.629(6) meV and Ezb2 = 
3.599(13) meV, respectively. The peaks observed in the 
constant wave-vector scans at Q^bi ^^^ Q2b2 ^^^ reso- 
lution limited. This experimental fact together with the 
identity of the values E^bi and Ezb2 confirms the NN 
interactions in 6c-plane are identical along the square- 
lattice directions. In case of different strengths for the 
NN interactions in 6c-plane a broadening of the peaks 
at Qzbi ^^d Q2b2 would be observed. From our pre- 
vious study^ we know that there is also a small NNN 
interaction equal to 2% of the NN exchange interaction. 
Therefore the observed one-magnon mode was compared 
to the following model 2D Hamiltonian: 



H 



(iJ> 



-jf(sr 



sr- 



-sf-sf)}+j2VSi-s 



i,k) 



(1) 

where (i,j) indicates the sum over NN in the 6c-plane, 
(i, k) - the sum over NNN in the 6c-plane, Jf, J^"^ and 
J2 are z-, x?/-components of the NN interaction and the 
NNN exchange, respectively. 
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FIG. 9: (Color online) The constant Q-scans collected at 
small energy transfer show the energy of a gapped spin- 
wave at the antiferromagnetic zone center performed at Q = 
(0,0,1) and at Q = (0,-1,0) are presented in (a) and (b), 
respectively, (c) Constant Q-scans performed close to the an- 
tiferromagnetic zone boundary at high energy transfer show 
the dispersion of the spin- wave. The measurements were per- 
formed in zero magnetic field and at T = 1.42 K. The solid 
lines represent the convolution of a Gaussian with the resolu- 
tion function. 



The linear spin wave theory (for details see Appendix 
B) yields two spin- wave modes with the dispersion hcjq = 



A^ 



B2 



where Aq = 4SJf + S(Jf - Jf )(cos(qb) 
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cos(qc)) — 4SJ2 + 4SJ2 cos(qb) • cos(qc)) and Bq = 
S(J^^ + Jf )(cos(qb) + cos(qc)). This imphes that the 
exchange anisotropy mostly affects the magnon energy 
close to the antiferromagnetic zone center, while the zone 
boundary energy remains nearly unaffected by the ex- 
change anisotropy. In the 2D S=l/2 AF, the energy 
of a classical (large- S) spin- wave mode is renormalized 
due to quantum fluctuations with the best theoretically 
predicted renormalization factor Zc = 1.18^^. There- 
fore the energy at the antiferromagnetic zone boundary 
is equal to Ezb = 2ZcJf — J2R, where J2R is the renormal- 
ized NNN interaction. The calculated xy-component of 
NN and NNN exchange interactions are equal to J^"^ = 
1.563(13) meV and J2 — 0.02J^^, respectively. Accord- 
ing to the linear spin wave theory E^^ = 8J^"^(J^"^ — Jf ) 
and thus Jf = 0.9979(2)Ji^. 

The values of the xy- and 2:-components of the NN in- 
teraction obtained from our neutron measurements are in 
a good agreement with the result of magnetic susceptibil- 
ity measurements^^, which yielded J^"^ = 1.507(26) meV 
and Jf = 0.9954J^"^. The small XY anisotropy indicates 
that the dominant exchange interaction between nearest 
copper ions in 6c-plane in Cu(pz)2 (0104)2, while spa- 
tially very anisotropic, is close to the isotropic limit in 
spin space, explaining the strong similarity of the HT 
phase diagrams measured in magnetic fields applied par- 
allel and perpendicular to copper square-lattice (Fig. [3]) . 

The inelastic-scattering data were fitted with the 
Oaussian instrumental resolution function convoluted nu- 
merically with the model Hamiltonian Q. The result of 
the fits is shown by the red lines in Fig. |8[and Fig. |9) and 
provides a good description of the observed spin waves. 
The color plot of the neutron scattering intensity, which 
is shown in Fig. [Toj summarizes the observed magnetic 
excitations in both crystallographic directions. The black 
lines display the result of the linear spin wave theory, 
showing that the observed dispersive excitation is well 
characterized by the Hamiltonian (IT]). 

The measured spin wave dispersion is similar to that 
observed in another 2D square-lattice antiferromagnetic 
material, namely copper deuteroformate tetradeut cr- 
ate (OFTD), where the exchange interaction strength 
is equal to J = 6.3(3) meV and the energy gap 
of E = 0.38 (2) me V is present at the antiferromag- 
netic zone center. However, the energy gap in OFTD 
is induced by the presence of small antisymmetric 
Dzyaloshins ky-Mo riya interaction D = 0.0051(5) meV 
between NN^^. Another example with comparable 
properties is K2V3O8 with 2D NN exchange strength J = 
1.08(3) meV and small energy gap at antiferromagnetic 
point equal to E = 0.072(9) meV^, which is described by 
Dzyaloshinsky-Moriya and easy-axis anisotropics^^. In 
contrast, Dzyaloshinsky-Moriya interactions between NN 
in Ou(pz)2 (0104)2 are forbidden by symmetry and the 
energy gap at the antiferromagnetic zone center is gen- 
erated by small XY anisotropy. 

We also studied the spin-wave dynamics along the an- 
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FIG. 10: (Color online) Color plot of the scattering inten- 
sity, showing the dispersion along the Q = (0,-1,/) and 
Q = (0, /c, 1) directions measured at zero field, presented 
in the left and right panels, respectively. The color plot 
was obtained by merging a total of five and thirteen con- 
stant wave-vector scans, respectively. The solid line repre- 
sents the dispersion computed from linear spin wave theory 
using Jf = 1.563 meV, J? = 0.9979Jf and NNN exchange 
equal to J2 = 0.02Jf as described in the text. 



tiferromagnetic zone boundary by performing constant 
wave-vector scans along the Q = (0,0.5,/) direction 
from / = 1 to / = 2. Typical data are shown in 
Fig. pT]^a-c) and the observed zone boundary dispersion 
is shown in Fig. pTJ^d). The onset of the scattering at 
Q = (0,0.5,1.5) is reduced by 10.7(4)% in energy com- 
pared to Q = (0,0.5,1), confirming our recent indepen- 
dent measuremenlP^. The decrease of the resonant mode 
energy at Q = (0,0.5, 1.5) results from a resonating va- 
lence bond quantum fluctuations between NN spins^^^. 
The observed dispersion at the zone boundary is slightly 
larger than expected from series expansion calculations 
and Quantum Monte Carlo simulations for 2D Heisen- 
berg square-lattice AF with NN interactions^ ^ and can 
be explained by the presence of a small antiferromagnetic 
NNN interaction. 

In order to subtract a nonmagnetic contribution from 
the background in the energy scan at Q = (0, 0.5, 1.5) we 
performed measurements with the sample turned away 
from magnetic scattering. The background-subtracted 
data are shown in Fig. pTj^a). The width of the scatter- 
ing peak as a function of energy at Q = (0,0.5,1.5) is 
clearly broader than the instrumental resolution. This 
implies the existence of a magnetic continuum scattering 
in this region of the antiferromagnetic zone boundary. 
The magnetic continuum with the present PANDA mea- 
surements is consistent with our previous investigation 
of Ou(pz)2 (0104)2^^. This non-trivial magnetic contin- 
uum and the dispersion at the zone-boundary result from 
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FIG. 11: (Color online) (a) The energy scan at Q = 
(0,0.5,1.5) with the background subtracted as explained in 
the text. The energy scans performed at wave- vectors Q = 
(0,0.5,1.7) and Q = (0,0.5,2) are shown in (b) and (c), re- 
spectively. The red curves are fits of a Gaussian function con- 
voluted with the resolution function, (d) The antiferromag- 
netic zone boundary dispersion measured in zero magnetic 
field and at T = 1.42 K. 



ture by neutron diffraction shows 2D antiparallel align- 
ment of spins at the neighbor sites. 

At zero field, we observed a well-defined magnon mode 
using neutron spectroscopy. Its dispersion is described 
by the spin Hamiltonian with slightly anisotropic NN 
interaction with xy- and z-components equal to J^^ = 
1.563(13) meV, J^ = 0.9979(2)Ji^ and NNN exchange 
equal to J2 — 0.02J^"^. Therefore the closeness of the 
dominant exchange interaction Ji to the isotropic limit 
measured by neutron spectroscopy explains the similar- 
ity of the HT phase diagrams obtained by specific heat 
measurements in magnetic fields applied parallel and per- 
pendicular to 6c-plane. 

Our results obtained by three different experimental 
techniques confirm and supplement each other clearly 
demonstrating that Cu(pz)2 (0104)2 is the first weakly 
frustrated 2D S=l/2 AF on a square lattice with the ab- 
sence of Dzyaloshinsky-Moriya interaction between NN. 
The measurements verify a relatively large 10.7(4)% 
zone-boundary dispersion and a rather strong magnetic 
continuum at the zone boundary. We associate these fea- 
tures with a resonating valence bond fluctuations which 
are enhanced by a small NNN AF interactions. 



quantum fluctuations in Ou(pz)2 (0104)2 which are en- 
hanced by a small antiferromagnetic NNN interaction. 



IV. CONCLUSIONS. 
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We have performed a comprehensive study of the 2D 
S=l/2 square-lattice AF Cu(pz)2 (0104)2 using a collec- 
tion of experimental techniques, including specific heat 
measurements, neutron diffraction and neutron spec- 
troscopy. The HT phase diagram was mapped out for 
magnetic fields up to /ioH = 9T applied parallel and 
perpendicular to the Cu^+ square-lattice planes, show- 
ing that the HT boundaries of the ordered phase have 
the same field dependence. This result shows that the 
dominant exchange interactions between nearest spins 
are close to the isotropic limit. Applied magnetic fields 
induces XY spin anisotropy leading to quench of quan- 
tum fluctuations, as we observed in the specific heat mea- 
surements. 

Neutron diffraction confirms the HT phase diagram 
obtained from macroscopic measurements and extends 
it up to /ioH = 13.5 T. The ordered magnetic mo- 
ment of Cu^+ ions at zero field, mo = 0.47(5) /iB, is 
reduced from the expected value, indicating the exis- 
tence of strong quantum fluctuations in the ground state. 
Magnetic fields quench quantum fluctuations and signifi- 
cantly increase the size of the ordered magnetic moment 
to mo = 0.93(5) /ie at iiqR = 13.5 T. Therefore neutron 
diffraction confirms the suppression of quantum fluctua- 
tions by field-induced anisotropy seen in our specific heat 
measurements. The study of the ordered magnetic struc- 



VI. 



APPENDIX A: THE GROUP THEORY 
ANALYSIS. 
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TABLE II: The character table and the irreducible repre- 
sentations obtained by performing group theory analysis for 
monoclinic space group C2/c (#15), the table setting choice 
is bl) and the magnetic ordering vectors k = (0, 0, 1) and 
k= (0,0,0). 

The crystal structure of Cu(pz)2 (0104)2 belongs to 
the monoclinic C2/c space group (#15), whose Laue 
class and the point group are 2/m. The Cu^+ ions 
occupy 4e Wyckoff positions and they are located at 
ri = (0 0.7499 0.25), r2 = (0 0.2501 0.75), T3 = 
(0.5 0.2499 0.25) and r4 = (0.5 0.7501 0.75). For 
C2/c space group, reciprocal lattice points are located at 
Q = (/i, /c, /) with h-\-k = 2n. The magnetic Bragg peaks 
were observed at Q^^ = (0,1,0) and at Q2 = (0,0,1) 
indicating that only Q2 can be the magnetic ordering 



vector. However, taking into account nearly identical 
crystallographic parameters b and c and the presence of 
be and cb crystallographic domains, the magnetic Bragg 
peaks at Q^ and Q2 are indistinguishable. Therefore we 
can not identify whether the magnetic ordering vector is 
k = (0, 0, 1) or k = (0, 0, 0). The subgroups of the order- 
ing wave vectors are identical and consist of four sym- 
metry operations which belong to four different classes: 
1, 25,1 and rriac' Here, 1 is the identity, 2^ is a two- fold 
rotation around the b-axis, 1 is the inversion and rriac is 
a mirror plane in the ac plane. Therefore, there are four 
one-dimensional irreducible representations whose char- 
acters are summarized in the character table given in 
Tab. [Hj The decomposition equation for the magnetic 
representation is r^ag = ITi + ir2 + 2Ts + '2T4. The six 



basis vectors presented in Tab. HI are calculated for two 
Cu^+ positions in primitive unit cell using the projection 
operator method acting on a trial vector (j)a 



\^^ 



E 

g€Gk 



Bt*{g)Y,Si,ei^s<f>adet{R^), 



where ^^^ is the basis vector projected from the A*^ row 
of the u^^ irreducible representation, D^*(g) is the A*^ 
row of the matrix representative of the u^^ irreducible 
representation for symmetry operation g, i denotes the 
atomic position and Rg is the rotational part of the sym- 
metry operation g. Note that basis vectors are identical 
for both k = (0, 0, 1) and k = (0, 0, 0). 

Dzyaloshinsky-Moriya (DM) interactions are defined 



as 



H 



DM 



= Y^ Dij-[SixSj], 



(iJ> 



where Dy is an axial vector. Action of any symmetry op- 
eration (including lattice translations) A on a DM vector 
Dij must be equal to DA(i)A(j) and Dy = — Dji. We an- 
alyze the action of the inversion symmetry operation 1 
on the axial DM vector D12, where i = 1 and j = 2 de- 
notes the NN copper positions ri = (0 0.7499 0.25) and 
r2 = (0 0.2501 0.75), respectively. The result of operation 
is 

T(Di2) = (D^2 D^:^ Bl,). 

The application of the inversion symmetry on the ions 
positions leads to 



1(D 



12 j 



D 



1(1)1(2) 



(Dfi DL D21 



These relations imply (Df^ D^2 Df2) = (D21 D^^ D|i) 
which is possible only in case of D12 = 0. Therefore, DM 
interactions between NN in Cu(pz)2 (0104)2 are forbid- 
den by the crystal symmetry. 



VII. APPENDIX B: THE LINEAR SPIN WAVE 
THEORY. 
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(0 10) 


(0 10) 
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(0 10) 
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Ta 


^3 
^4 


(10 0) 
(0 1) 


(10 0) 
(0 1) 


r4 




(10 0) 
(0 1) 


(-1 0) 
(0 -1) 



TABLE UL Six basis vectors calculated for two Cu^^ posi- 
tions in primitive unit cell as explained in Appendix A. 



tion we can make the Holstein-Primakoff transformation 
of the spin compounds into bosonic creation and anni- 
hilation operators. In linear spin wave approximation it 
gives: 

0^ =0 a^ a^, b; = o -p a- aj, bj^ ^ b aj^ aj^. 



Sf = V2S ^^^, S 



y 



V2S^L^ 



, — ai -\- a- , — ai -\- a- 



^k + \ ny _ ^A^ ^k - a,^ 

^2~' ^k-V2b^^^. 

The quantization axis lies in 6c-plane and therefore Jf = 
Jf = J, Ji = J — A. The spin Hamiltonian written in 
bosonic operators is: 



H = ^{4JS(a+ai + a+aj) + 2S ^^^^ (aiaj + a+a+)- 

(iJ> 



2S-(aia++a+aj)}+ 



y^ 4J2S{(aia^ + a+ak) - (a+ai + a^ak)}, 

where (i,j) indicates the sum over NN in the bc-plane, 
(i, k) - the sum over NNN in the bc-plane. After Fourier 
transformation obtained Hamiltonian can be diagonal- 
ized using standard Bogoliubov transformation: 

aq = -Uqaq+Vq/3+, 



aJ" = -UqOfJ" + Vq/3q 



Assuming the system in the antiferromagnetic Neel 
ground state with spins pointing along z and —z direc- 



^— q — ^q^q UqPq, 



10 



a-q = Vqaq-Uq/3+, 

where Ofq and /3q are the bosonic operators and Vq, Uq 
are numbers. Finahy, after the diagonahzation we have 

Hq = Eg.s. + ^S[2Aqa+aq + Bq(a+a+q + aqa_q)], (2) 



4SJ2 + 4SJ2 • cos(qb) • cos(qc) 



and 



and the eigenstates are given by hcjq = (A^ — Bq)-*^/^, 



where 



Bq = 2S(J - l/2A)(cos(qb) + cos(qc)). 



Aq = 4SJ + SA(cos(qb) + cos(qc))- 
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